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Analytical Solutions of Basic Models in 
Quantum Optics 


Daniel Braak 


Abstract The recent progress in the analytical solution of models invented to de¬ 
scribe theoretically the interaction of matter with light on an atomic scale is re¬ 
viewed. The methods employ the classical theory of linear differential equations in 
the complex domain (Fuchsian equations). The linking concept is provided by the 
Bargmann Hilbert space of analytic functions, which is isomorphic to the 

standard Hilbert space for a single continuous degree of freedom in quantum me¬ 
chanics. I give the solution of the quantum Rabi model in some detail and sketch 
the solution of its generalization, the asymmetric Dicke model. Characteristic prop¬ 
erties of the respective spectra are derived directly from the singularity structure of 
the corresponding system of differential equations. 

Key words: Quantum optics, Bargmann space. Differential equations. Singularity 
theory, Integrable systems 


1 Introduction 


The interaction of matter with light forms a major subject of theoretical and ap¬ 
plied physics HI. It is essentially characterized by the quantum nature of both con¬ 
stituents, studied within Quantum Optics ||2|. The quantum features of the processes 
occurring in recently realized nano-sized devices can be used to control the gener¬ 
ation of entangled states 0, thereby allowing to construct the basic elements of 
quantum information technology IH. 

The technological advances in nanofabrication made it possible to reach very 
large coupling strengths between the light (usually confined to a single or few modes 
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in a cavity) and an (artificial) atom described by a discrete set of energy levels in 
the deep quantum limit Ellllll. The prototypical system consists of a “matter” part 
with two possible states coupled to the dipole component of a single radiation mode. 
The Hamiltonian of the atom can thus be expressed via Pauli spin matrices O) and 
the radiation through a harmonic oscillator with frequency co. The Hilbert space of 
the total system reads then (8)L^(K) and the Hamiltonian 

Hu = (Oa^a +gGxia + a^)+ ( 1 ) 

Here o' and a are the creation and annihilation operators of the bosonic mode and 
energy is measured in units of frequency (h = 1). cOq denotes the energy splitting 
of the two-level system (the “qubit”) which is coupled linearly to the electric field 
(^ {a + a')) with interaction strength g. This model was studied semiclassically 
already in 1936 by Rabi [Sl and the fully quantized version Q has been introduced 
in 1963 by Jaynes and Cummings ||9l. It is therefore called the quantum Rabi model 
(QRM). Despite its apparent simplicity, the QRM is difficult to solve analytically 
because it does not exhibit invariant subspaces of finite dimension like the following 
model, 

Hjc = coa^a +g{a^a +G^a^)+ ( 2 ) 

with = ((7;c±i<Tv)/2, which corresponds to the “rotating-wave” approximation 
of 0 0. The Jaynes-Cummings model (JCM) 0 can be justified close to reso¬ 
nance, 0) ~ Wo and small coupling g/co < 1 Go] and has been the standard model 
for typical quantum optical applications with g/(0 < 10^^ for many years. The ma¬ 
jor simplification arising in 0 as compared to 0 consists in the fact that the op¬ 
erator C = a^a -b commutes with Hjc, which means that decays into 

infinitely many Hyc-invariant subspaces in which C takes constant values G INo, the 
set of non-negative integers. Each of these spaces is two-dimensional entailing triv¬ 
ial diagonalization of Hjc- The fact that the polynomial algebra of C has infinite 
dimension means that C generates a continuous t/(l)-symmetry of Hjc ifTOl fTTIl : 
defining t/(0) = exp(/0C), we have {(j))aU{(p) = e“l’a, U\(j))a'^U{(j)) = 
and W{<p)a^U{<p) = This abelian symmetry associated with the integral 

of motion C renders the JCM integrable, because it has only two degrees of freedom, 
the continuous one of the radiation mode with Hilbert space L^(M) and the discrete 
one of the qubit with Hilbert space C^. 

The concept of integrability underlying this argument amounts to a direct transfer 
of Liouville’s definition from classical mechanics to quantum mechanics: A system 
with N degrees of freedom is integrable if it exhibits N independent phase space 
functions which are in involution with respect to the Poisson bracket lfT2l . In the 
JCM these are the Hamiltonian Hjc and C. But because independence of opera¬ 
tors cannot be defined in analogy to functions on phase space, this definition is not 
feasible as any Hamiltonian system would be integrable according to it Ha. 

The continuous t/(l)-symmetry of the JCM is broken down to a discrete symme¬ 
try by the counter-rotating term (7+in the QRM. U^7z)HrU {n) —Hr and 
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Hr commutes with P = (—1)" “(7j = —U{%). Because = 1, its polynomial alge¬ 
bra is two-dimensional and P generates a Z 2 -symmetry of Hr, usually called parity. 
The eigenvalues ± 1 of P characterize two ///^-invariant subspaces (parity-chains), 
each of them infinite-dimensional m. Therefore, the problem appears to be only 
marginally simplified by using the parity symmetry and it was widely held that the 
QRM is not integrable im. However, it could be demonstrated that the weak par¬ 
ity symmetry is indeed sufficient for integrability of the QRM because it possesses 
only one continuous degree of freedom, whereas the Hilbert space dimension of the 
discrete degree of freedom matches the dimension of the polynomial algebra gener¬ 
ated by P, rendering the QRM integrable according to the level-labeling criterion for 
quantum integrability El . The detailed understanding of the spectrum (and dynam¬ 
ics) of the QRM beyond the rotating-wave approximation has been necessitated by 
the recent experimental access to the ultra-strong and deep-strong coupling regime 
in circuit QED 0121 and through quantum simulations. 

This survey is organized as follows: In section]^ the analytical solution of the 
QRM is presented based on a formulation of the problem in Bargmann’s space of 
analytical functions, section deals with multi-qubit models and the last section 
contains some remarks on possible future research directions. 


2 The Quantum Rabi Model 

The Z 2 -symmetry of the QRM can be used to eliminate the discrete degree of free¬ 
dom from the problem, just as the f/(l)-symmetry of the JCM allows elimination 
of the continuous degree of freedom. Each parity-chain is isomorphic to L^(IR) 
and Hr reads in 


+ (3) 

with A = 0)0/2. The complication of this reduced Hamiltonian comes from the 
last term (—1)" On the other hand, this term is instrumental for the analytical 
solution of the model. To elucidate its meaning, it is convenient to represent (|^ 
in Bargmann’s space of analytical functions which is isometrically isomorphic to 
l2(k) ITb). The space SS is spanned by functions /(z) of a complex variable z 
which have finite norm (0|0) with respect to the scalar product 

{\l/\(j)) = - [dzdze~^\i/{z)(j){z) (4) 

n J 

and are analytic in all C (dzdz = d9I(z)d3 (z)). The criterion for being an element of 
£§ is therefore two-fold: (j) (z) G ^ if both of the following conditions are satisfied: 
(B-I): (010) <oo 

(B-II): 0(z) is holomorphic everywhere in the open domain C. 

The isometry maps f{q) G L^(R) to an analytic function 0(z) € 
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Hz) = Af]iz) = ^ fjq (5) 

The operators H,a are mapped to z and d/dz, respectively, 

= (6) 

dz 

The normalized vacuum |0) with a|0) = 0 is mapped to the constant function Hiz) = 
1. On infers from Q that all polynomials in z are elements of Especially the 
n-th eigenstate of the harmonic oscillator \n) ^ I^Hn{q) is mapped onto the 
monomial f' j^fn}.. Moreover, all functions which have the asymptotic expansion 

0(z) =e“‘^z^““(co+ciz^'+C 2 Z^^ + ...) for z-)-°o, (7) 

with arbitrary ai € C satisfy (B-I) iflTl . (0 is the asymptotic form of the normal 
solutions of a differential equation having an unramified irregular singular point of 
s-rank two at infinity Haim. The coherent state |a) = |0) in L^(IR) is mapped 

obviously to the exponential function Functions behaving as 

(j){z) =exp|^^z^ + aiz^z^““(co + ciz^' +C 2 Z^^ + ...) (8) 

asymptotically, satisfy (B-I) only if |o! 2 | < 1 and correspond to an irregular singu¬ 
larity of s-rank three. The limiting value \a 2 \ = 1 belongs to the (not normalizable) 
plane wave states fp{q) — exp(/pq')/ Hln, 

■^[fp\iz) = ^^^exp [^z^ + iVlp^ . (9) 

The Hamiltonian H+ reads in (with 0 = 1), 

H+=z^z+s[z+^^+^f, (10) 

where T denotes the reflection operator r[v/](z) = y/(—z). The Schrodinger equa¬ 
tion (//+ — E)\j/{z) =0 corresponds to a linear but non-local differential equation in 
the complex domain, 

z^r(z)-fg(|^-fz) \i/{z)=Ey{z)-Ay{-z). (11) 

The theory of these equations initiated by Riemann and Fuchs ifTSll can now be 
applied to First, one obtains with the definition v/(z) = 0i(z) and i/(—z) = 
H (z) the coupled local system. 
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iz + g)^<l)iiz) + igz-E)(l)i{z)+A(l) 2 {z) =0, (12a) 

iz-g)^(j)2{z)-igz + E)(j)2iz)+A(j)i{z) = 0. (12b) 

az 

This system has two regular singular points at z = ±g and an (unramified) irregular 
singular point of s-rank two at z = oo QSl. The normal solutions of ( fT^ behave 
asymptotically as Q with ai = g or ai = —g. The two Stokes rays are the positive 
and negative real axis. We infer that all solutions of ( |T2l i satisfy (B-I). I follows 
that the discrete spectrum {£„}, n G INq will be determined by (B-II), because not 
all solutions of ( [I^ are analytic in C. Define x = £ + g^. Then the exponents of 
(pi (z) at the regular singular point g (—g) are {0, 1 +x} ({0,x}), while for (p 2 {z) the 
exponents atg (—g) are {0,x} ({0,1 +x}). 

If £ belongs to the spectrum of (pi and (p 2 must be analytic in C, especially 
at both points ±g. This leads naturally to a division of the spectrum into two parts: 

1) The regular spectrum o>e^ consisting of those values £„ for which x„ = £„+ g^ 
is not a non-negative integer. 

2) The exceptional spectrum Gexc for which Xn G INq. 


2.1 The Regular Spectrum 


\i x^ INq, the only allowed exponent at both points ±g is 0. We consider the case 
—g and define y = z + g, (pi .2 = 0i.2. Then, 


y—(pi = x(pi-A(p 2 , 
dy 

(y-2g) — 02 = {x-4g^+2gy)(p2-A(pi. 
dy 


(13) 

(14) 


A local Frobenius solution for (p 2 (y), analytic at y = 0, reads 02 (y) = T,7=o^n(^)y'‘ 
with coefficients K„(x) to be determined. Integration of ( |T3] l yields 

^i(y) =c/-^f^K„(x)^^. (15) 

«to 


Because x ^ INq, c must be zero. This determines 0i (z) uniquely in terms of 02(z). 
Setting Ab = 1, the following three-term recurrence relation for the K„ (x) is obtained 
from ( firl i, 

nKn = fn—l{x)Kn—i Kn—2t ( 16 ) 

with 

/„(X) =2g-b;^ 



x — n 


(17) 




6 


Daniel Braak 


and initial condition Kq — l,/ri(x) = ^liy) will usually develop a branch- 

cut singularity at y = 2g. The radius of convergence of the Frobenius solution 
around y = Qis R = 2g, which can be deduced from the asymptotic value 1 / (2g) of 
fn-\{x)/n for n —>■ oo. Due to the relation (j) 2 (z) = V^(—z), the formal solution of ( |T3] l 
and ( [l4| i yields two expansions for V^(z), one analytic at z = g and the other analytic 
at z = -g, 


oo 


V^(z) = 02(-z) = ^ K„{x){-z + gY, 

(18) 

Wiz) = 01 (z) = £ K„ {x)A + 

(19) 


«=0 


The two circles centered at z = ±g are shown in Fig. Because the vectors 


Fig. 1 The singularity struc¬ 
ture of {n} and B21 . Two 
local Frobenius solutions an¬ 
alytic at z = g and z = —g 
respectively are defined by 
(Tt). If they coincide 
in the intersection of their 
domains of convergence, they 
describe the same (analytic) 
function in C. 



((pi(z),(p 2 (z))'^ and ( 02 (—-z), (—satisfy both the homogeneous hrst-order sys¬ 

tem they coincide in a neighborhood of zo if 

01 (zo) = 02(-Zo), 01 (-zo) = 02(zo) (20) 


for any zo in the intersection of their domain of convergence. That means that 02 (—z) 
is the analytic continuation of 0i(z) and itself analytic at z = g, therefore ij/(z) is 
analytic at both singular points. Both conditions in ( |20| ) are equivalent if zo = 0 
1201. This leads to the dehnition of the G-function for the regular spectrum of 

El EH, 


G+(x) = M0)-^i(0)=’£k„(x) 1 

fi=0 



( 21 ) 


If G+[En g^) = 0, the corresponding formal solution analytic every¬ 

where and an element of because it satishes (B-I) and (B-II), entailing that 
En G Greg{H+). G_(x) for H is obtained from G+(x) by replacing A with —A 
in ( |2T] i. It follows from ( [T6| ), ( [T7] i and (|2T| that G±(x) has simple poles at x G INq. 
The zeros of G±(x) are distributed between these poles. Fig. shows G±(x) and 
Fig. [^the corresponding spectrum of Hr for both parities. It is easy to see that the 
regular spectrum is never degenerate, neither within each parity chain nor among 
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Fig. 2 The G-functions for 
odd (blue) and even (red) 
parity for g = co = 1 and 
4 =0.4. 



states with different parity m . The G-function can be expressed in terms of known 
special functions as follows 

G±(x)= r,S,p,a-,1/2)a-,1/2). (22) 

Hc{a,Y,5,p,G',z) denotes a confluent Heun-function fT2\ and //'(a, 7 ,5,p,(J;z) 
its derivative with respect to z- The parameters are given as 

a =—X, 7=1— .r, d = —x, 

P = -g^, a=x{4g^-x)+A^. 

The functional form ( |2T] l of G± (x) leads to the following conjecture about the dis¬ 
tribution of its zeros along the positive real axis. 



Fig. 3 Rabi spectrum for the same 4 , (U as in Fig.|^and 0 < g < 0.8. The intersections between the 
spectra of different parity indicate the degenerate part of the exceptional spectrum. The two-fold 
labeling of states on the left corresponds to the uncoupled system (the ± denotes the spin quantum 
number) and on the right to the coupled case (± denotes the parity quantum number). 
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Conjecture; The number of zeros in each interval [n,n + 1], n G INq is restricted to 
be 0, 1, or 2. Moreover, an interval [n,n + 1] with two roots of G±{x) = 0 can only 
be adjacent to an interval with one or zero roots; in the same way, an empty interval 
can never be adjacent to another empty interval. 


2.2 The Exceptional Spectrum 

We shall demonstrate in the following that the presence of the exceptional spectrum 
Oexc poses certain constraints on the model parameters g and A such that for given g, 
A, at most two eigenvalues are exceptional. Furthermore, Gexc = (^exc C where 
comprises the values E = m — g^ with m G IN. Each eigenvalue in is doubly 
degenerate among states with different parity, is not degenerate and may take 
values E = m — g^ with m G INq. 

We begin with d"^. The poles of G+(x) at x S INq indicate that an integer x can 
only signify an eigenvalue of //+ if the corresponding pole in G+(x) is lifted for 
special values of the parameters g and A. If x € INq, not only the exponent 0 but 
also the exponents x, respectively x + 1 guarantee analyticity of </>i{z) and 02 (z) 31 
z = —g. However, as the difference of the two exponents at both singular points is a 
positive integer if x > 0 (for x = 0 this difference is positive at one singular point), 
the local analytic Frobenius solutions around z = —g will develop a logarithmic 
branch-cut at z = g in general. For x = m G INq, there exist always a solution for 
02 (y) analytic at y = 0 of the form 

Uy)= L Kny\ (23) 

n=m+\ 

because the largest exponent of 02 (y) aty = 0 is x-f 1 ifTSll . Integration of ( [T3] l yields 
for 01 (y), 

0i(y)=cy'"-4 £ K„^. (24) 

^,1 n — m 

n=m+L 

In this case, the constant c may be different from zero because 0i (y) is then ana¬ 
lytic at y = 0. Solving now with the ansatz ( | 2 ^ , we obtain for n>m + 2 the 
recurrence ( [Thl l, ( [T7| i and the initial conditions 

= K^ = 0. (25) 

2[m+l)g 

c is fixed in terms of ■ Setting = 1 we obtain for \j/{z) the two expressions 
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= 02(-z) = Y, K„{m;g,A){-z + g)", 

n=m+l 

V,(z) = + £ K,(,n-g,A)^-^±^ 

and the G-function follows as 


GW(,,a) = -?l 2 L^ 


- Y Knim;g,A)[l- 

n=m-\-\ 


^n—m— 1 


(26) 

(27) 


(28) 


The zeros of the function G^^ (g, A ) determine those values of the parameters g and 
A for which //+ has the exceptional eigenvalue m — g^ with m G INq. For odd parity, 
we have G^!"\g,A) = G^_^\g,—A). It follows that G^^^\g,A) and G^"'\g,A) have 
no common zeros, so this part of the exceptional spectrum is not degenerate, just as 
the regular spectrum. It was computed by a related method in ||23|. 

To obtain we consider now the smaller exponent, zero, of ^{y) at y = 0, 
leading to the expansion 

Uy)^YKny\ (29) 

n =0 

After integration of ([T3ll, (y) reads 


(j)l (y) = cy"’ -A Yk„ — - Ay”'Km ln(y). 

^ n — m 

n^m 


(30) 


The K„ for n < m — I are determined again with and initial conditions Kq = 1, 
Ki = fo{m). Therefore Km{m-,g,A) is uniquely hxed. The logarithmic term in ( [30l l 
vanishes if Km(ni-,g,A )=0ll5l. The coefficients Kn for n > m + 1 are determined 
with ([T6|) and initial conditions 


^m+l — 


m+ 1 


— -K 

J^m— 1 

^8 


Km = 0 . 


In this case there exist two local solutions analytic at y = 0, ( |2^ and ( |29] l. If m ^ 0, 
they are linearly independent and span the whole solution space for 02 (y)- Because 
of the reflection symmetry mapping y to 2g — y, these solutions describe the solution 
space in a neighborhood of y = 2g as well and all solutions of ( [T^ are analytic at 
both g and —g, thus in all of C, if Km(m',g,A) = 0 and no further condition is 
necessary. Moreover, the Km's are the same for odd parity, so the eigenvalue E = 
m — g^ obtained via this condition is always doubly degenerate between states of 
different parity. The presence of the spectrum and its “quasi-exact” nature m 
may be explained more generally in terms of the representation theory of s( 2 (K) 
1251 . A special situation arises forx = 0. The condition Kq{0) — 0 renders the ansatz 
( |29l l equivalent to ( |23| l and both solutions are linearly dependent. Thus only one 
local solution may be analytic at y = 0 and 0i (y) is given by (|24ll for m — 0. If 
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E = —g^ is an eigenvalue of H+, G^\g,A) must vanish. It follows that for these 
parameter values G^_^\g,A ) ^ 0, the eigenvalue E = —g^ is never degenerate and 
an element of cr'^. 


2.3 Methods based on Continued Fractions 

The Bargmann space formalism has been applied to the QRM as early as 1967 
by Schweber j^ . He did not make use of the Z 2 -symmetry but derived the cou¬ 
pled system ( [T^ directly from Q with the ansatz ^ = {(pi(z),(p 2 iz))^ € 
for the wave function with energy E. ( [T^ is then satished by 01 = + (p 2 and 

02 = ^1 — <P 2 - He obtained the local Frobenius solution for 02 (y) given by ( [Thl l and 
The convergence radius of the series ( |29] l is 2g for arbitrary x and the discrete 
set of eigenvalues is selected by determining those x for which ( |29] l has inhnite con¬ 
vergence radius. The problem is equivalent to compute the minimal solution of the 
recurrence ( [Thl l ll27ll . The spectral condition obtains then by equating the minimal 
Ki {x)/Kq with fQ{x) from the initial conditions. The equation has the form T’(x) = 0, 
where T'(x) is represented by a continued fraction ESI . This method, while formally 
correct, has several conceptual shortcomings: 

• The function F{x) has an unknown singularity structure and it is impossible to 
infer qualitative aspects on the distribution of its zeros from it. 

• The actual computation of the continued fraction makes a truncation at some 
order necessary which is equivalent to define the model on a finite-dimensional 
Hilbert space, which is the starting point for other work employing continued 
fractions to compute the Rabi spectrum 1^12^ . 

• The zeros of F{x) correspond to Greg U but there is no possibility to discern 
both types of spectra, especially the double degeneracy of cannot be detected 
with this method. 

• The spectrum G"fg is not accessible because the expansion ( |29l l with Ab 7 ^ 0 is 
assumed in the derivation of F{x). 

Besides these conceptual problems, the method is numerically feasible only for 
the hrst low-lying eigenvalues. The continued fraction has a pole in close vicinity 
of each zero and their distance approaches zero exponentially for growing x, so that 
at most ten energy levels may be resolved within a double precision calculation. On 
the other hand, the equivalence of the continued fraction approach to exact diago- 
nalization in finite-dimensional Hilbert spaces proves the validity of the latter for 
the QRM flTl . 

Schweber’s technique is confined to problems reducible to three-term recurrence 
relations for the local Frobenius solutions and implements then (B-II) as the spectral 
condition. It fails for models with more than a single qubit because the ensuing 
recurrence relations have more than three terms. The next section is devoted to the 
application of the theory presented above to models with > 1 qubits. 
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3 The Dicke Models 


The natural generalization of the Hamiltonian Q couples several two-level systems 
to the same mode of the radiation field, 

A ^ N 

Hdn = COa^a + V -7r<yiz + + ( 31 ) 

ti 2 VNti 


This model assumes different qubit frequencies cuo; and couplings gi to the field and 
is therefore called the asymmetric Dicke model (ADMat) with N qubits. Dicke in¬ 
troduced the permutation-invariant version of iD in 1954 and predicted the (later 
observed) phenomenon of “superradiance” for large N 1^ . Its rotating-wave ap¬ 
proximation is integrable for all N ED, while the full model is non-integrable for 
any A > 1 according to the level-labeling criterion ifTSll . Applications in quantum 
information technology mandate the study of ( [STj l without approximations for small 
N, because it describes the implementation of quantum gates within circuit QED 

ED- 


The following section treats the asymmetric model for N = 2 and section 3.2 the 
symmetric model for A = 3. 


3.1 ADM 2 and Exceptional States 

The Hamiltonian of the ADM 2 reads in slightly different notation, 

Hd2 = (Oa^a -\-gi(7ix{a-\-ci') g2(^2xio- -f AiCij + A 2 G 2 z- (32) 

This model has a Z 2 -symmetry similar to 0 ’ generated by P = exp{i7ta‘a)(Ji^(J 2 z- 
However, because it has only two irreducible one-dimensional representations, the 
discrete degrees of freedom cannot be labeled with a Z 2 -quantum number (their 
Hilbert space is C'^). The symmetry of the model is not sufficient to make it inte¬ 
grable as the A = 1 case 0 . 

Nevertheless, the same methods as above can be used to solve ( |32| l exactly ll33ll . 
at the expense of a more complicated G-function, which is no longer a linear combi¬ 
nation of formal solutions as ( [2T| . After application of the symmetry, the remaining 
Hilbert space is not but The Hamiltonian reads (o) = 1), 

+ (^1 +82(yz) (^+ ^) + {^2±Ait)ax, (33) 

with f = (—1)^32. An eigenfunction of ( |33| l with eigenvalue E is the vector yf = 
((pi {z),(p 2 {z)y. Defining ^ 3 ( 2 ) = (pi (—z) and (p 4 {z) = (pzi—z), we obtain a coupled 
system of four ordinary first order differential equations with four regular singular 
points located at g = gi +g 2 , g’ = gi— g 2 , —g and —g'. Moreover, z = 00 is an un- 



12 


Daniel Braak 


ramified irregular singular point of s-rank two; we conclude that again all formal 
solutions fulfill (B-I) and E is determined by postulating analyticity of the solution 
at all regular singular points. The corresponding coupled recurrence relations for 
the Frobenius expansions around each of the points 0, g and g' cannot be reduced 
to a three-term recurrence, except in the case g' = 0, which allows a treatment sim¬ 
ilar to the QRM. It turns out that eight initial conditions determine functions (z), 
k = 1,...,32, describing the (Pj(z) around different expansion points. Overall an¬ 
alyticity is then equivalent to the vanishing of the determinant of a 8 x 8-matrix 
M±{E), whose entries are composed of the (j)/^, evaluated at the ordinary points zo 
and Zq, whose location depends on the geometry of the analytic regions 133]. The 
G-function can then be defined as G±{E) — det(M±(£’)). This function has poles at 
integer values of E +g and E +g , defining the exceptional spectrum, besides the 
regular, given by the condition G{E„) = 0. 

The spectra of Hj )2 obtained in this way are depicted in Fig. as function of g 
and various levels of asymmetry. For the completely asymmetric cases in Figs.Qa), 
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Fig. 4 The spectra of Hm with (a) Ai = 0.6, A 2 = 0.2, 0 ) = l,0<g = gi+g 2 < 2.5, gi = 4g2. (b) 
Ai = 0.6, A 2 = 0.2, CO = 1, 0 < g = gi -|-g 2 < 2.5, gi = 2g2. (c)4i = ZI 2 = 0.5, to = 1, 0 < g = 
gl + 82 < 2.5, gi = g 2 . (d)4i = 0.6, A 2 = 0.4, CO = 1, 0 < g < 2.5, gi = g 2 . Blue lines are 
eigenvalues with odd parity, while red lines are eigenvalues with even parity. 


0 b), we observe level crossings between states of different parity, whereas states 
with equal parity show avoided crossings, some of them quite narrow as the insets 
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demonstrate. There is no relation between degeneracies and the exceptional spec¬ 
trum as in the case of the QRM, because the G-fanctions G+ {E ) and G {E ) are 
not simply related. The set of lines with E + = n, n G INq, give the asymptotic 

(parity degenerate) spectrum in the deep strong coupling limit. FigQc) shows the 
completely symmetric case, Ai = A 2 , gi = g 2 - The invariance under permutation 
symmetry of ( [3^ leads to separation of the Hilbert space of the spin-1/2 qubits into 
singlet and triplet sector according to \ = 0©1. The total Hilbert space be¬ 
comes thus The singlet subspace is isomorphic to and 

the Hamiltonian describes the decoupled radiation mode. Thus the eigenenergies are 
just integer multiples of (£>, independent of the coupling. They are seen as horizon¬ 
tal lines in Fig. Qc). The triplet subspace is coupled to the radiation field and the 
spectrum shows a nontrivial dependence on g. 

An interesting situation obtains for equal couplings gi = g 2 but different qubit 
energies, Ai 7 ^ A 2 . The full permutation symmetry is broken, but there is a certain 
remnant of it. For gi = g 2 , z = 0 is a regular singular point and there exist “quasi- 
exact” eigenstates, belonging to the exceptional spectrum with £„ S IN for certain 
parameter values of g, Aj. These states contain a finite number N of photons (con¬ 
trary to the likewise quasi-exact elements of C 7 ^(, in the QRM) and are determined 
by a polynomial equation for g,Ai,A 2 depending on the energy value E = N. It 
reads for N = 2, 


2 - 


{A2±A\f 


{1-{A2T^i?)-8^ 


(tAi -A 2 ) =0, 


(34) 


where the + (—) corresponds to even (odd) parity. The condition comprises the 
symmetric case Ai = A 2 for odd parity, but otherwise determines Ai, A 2 in terms of 
g. This is true for all N >2. However, for A = 1 we find. 


(±Ai-A2)[1-(A2±Ai)2] =0. (35) 

There exists an eigenstate with energy E = I and even (odd) parity, if Ai -f A 2 = 1, 
(|Ai — A 2 I = 1). This state contains at most one photon and the condition for its exis¬ 
tence does not depend on g. It was first discovered by Chilingaryan and Rodrfguez- 
Lara ll34ll . With the notation ji//) = |n,si,S 2 ) for a basis element in where we 
have used the occupation number basis for the boson mode and sj G {g,e} denotes 
the state of the y'-th qubit, the exceptional state with even parity (Ai -f A 2 = 1 ) reads, 

\Ve) = ^ 2(Ai^ - |l,e,g) -f |l,g,e)^ , (36) 

with a normalization factor -yV. This state becomes the singlet state 

-^\l)®{\g,e)-\e,g)) (37) 

in the symmetric case Ai = A 2 . The spectrum for Ai + A 2 = 1 is shown in FigQd), 
exhibiting the quasi-exact state with £ = 1 as a g-independent line. In contrast to the 
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fully decoupled singlet state the state ( [36l l is strongly coupled to the radiation 
held, as its components depend on g. It is quite remarkable that states with hnite 
maximal photon number exist for arbitrary strong coupling without making the ro¬ 
tating wave approximation. This feature cannot be realized in the QRM, where each 
eigenstate contains always an inhnite number of photons. Due to its very simple 
structure, the state ( |36| ) could be useful for quantum computing applications, espe¬ 
cially as the condition for its existence depends only on the (easily controllable) 
qubit energies and A 2 and not on the coupling strength. Similar states are ex¬ 
pected to exist in all models ADMat with g, = g and even N. 


3.2 ADM 3 


The Hilbert space of the symmetric model ADM 3 with = g' = V^g, OJoi = COq = 
2 A may be splitted according toj(g)j(g )5 = j 0 j 0 |. Each of the two spin-1/2 
components are equivalent to the QRM, while the last component describes a single 
spin-3/2 coupled to the radiation mode with Hamiltonian, 

Hd 3 = a^a + 2Af^ + 2g{a + a^)Jx, (38) 

where and Jx are generators of SU{2) in the (four-dimensional) spin-3/2 repre¬ 
sentation. The Z 2 -generator has here the form P = exp{ina^a)R. R is an involution 
acting in spin-space as RJ^R = Jz, PJxP = — Jx- Application of P gives the following 
differential operator in each parity subspace. 


H±=z-+A 

dz 


0 V3 \ _ 
V3 ±2f ) ^ 



(39) 


Employing now the same machinery as in the previous section, we obtain again 
four coupled first order equations having regular singular points at ±g, ±3g and 
an irregular singular point (s-rank two) at inhnity ll35]l . Because of the stronger 
symmetry of ( |38] l compared to ( [3^ , the matrix M{E)±, whose determinant gives 
the G-function is only 6 x 6 . It contains 24 functions evaluated at the points zq = 
2g and Zq = 0. The poles of G±{E) (giving cj"//) are located at £ S IN and 
E 0 9g^ € INq. The curves determined by the latter set are also the limiting values 
for the spectrum for very large coupling. Eig. shows the spectral graph of the 
model as function of the coupling g. As in the N = 2 case, the degeneracies occur 
within the regular spectrum between states of different parity. However, the fact 
that we have now a determinant as G-function means that in principle degeneracies 
within the same parity chain are not excluded as the corresponding matrix M± [E) 
could have a higher-dimensional kernel at a specific energy. Up to now there is no 
numerical evidence for this scenario and the existence of these novel degeneracies 
is an open question. 


Analytical Solutions of Basic Models in Quantum Optics 


15 



Fig. 5 The spectrum of the Dicke model for even (red) and odd (blue) parity at A = 0.7 and for 
varying g. The y-axis shows x = E + g^. The ground state has odd parity as in the QRM. The 
two ladders of eigenvalues with different parity intersect within the regular spectrum. There are 
no degeneracies (but narrow avoided crossings) for fixed parity in this parameter window. Dashed 
lines denote the set E + 9g^ 6 INq and emerge as limiting values in the deep strong coupling regime 
1 - 


4 Conclusions 

We have seen in the previous sections that the classical theory of linear differential 
equations in the complex domain may be used to solve exactly elementary but im¬ 
portant problems in the field of theoretical and applied physics, contributing in this 
way to a better understanding of the basic models of circuit QED, which has been 
envisioned as a promising environment to implement devices capable of performing 
quantum computations. 

The mathematical technique relies on the Bargmann Hilbert space which 
allows to represent the Hamiltonians from Quantum Optics as differential oper¬ 
ators acting on functions of a complex variable z. Of central importance is here 
Bargmann’s two-fold spectral condition, which demands not only normalizability 
with respect to the scalar product but also analyticity in C for any function (j) (z) 
being an element of In this way it becomes possible to use the easily accessi¬ 
ble singularity structure of the corresponding differential equations to implement 
the spectral condition without recourse to a polynomial ansatz for the wave func¬ 
tions, which works for elementary integrable systems like the harmonic oscillator, 
the hydrogen atom and the Jaynes-Cummings model, but fails already for the quan¬ 
tum Rabi model, which is nevertheless integrable in a well-defined sense JBll (for a 
recent comparison with Yang-Baxter integrability see mi 

The method can be extended to models of central relevance for quantum technol¬ 
ogy, the Dicke models with a small number of qubits. These models are no longer 
integrable in view of the level labeling criterion iia but exactly solvable with the 
presented technique, although many characteristic simplifications of the quantum 
Rabi model are absent. Further applications concern models with a single irregular 
singular point as the two-photon Rabi model El, or the anharmonic oscillator. One 
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may also try to extend the formalism to multi-mode models ll38l . Notwithstanding 
these generalizations, we note that the already solved systems give rise to a set of 
open mathematical problems like the conjecture on the level distribution presented 
in section |2.1[ or the question whether a novel class of degeneracies exists in the 
Dicke models (section|3^. Thus, a future research direction will be the exploration 
of the recently observed connection 091 between the quantum Rabi model and the 
non-commutative harmonic oscillator Ho). It should be possible to transfer mathe¬ 
matical methods used for the study of the latter to the problems mentioned above. 
The techniques used in BTl 1^ l43l l44l could be applied f.e. to investigate the level 
crossing appearing in the ground state of the anisotropic quantum Rabi model ll45l . 
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